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WATER  DROPLET  EVAPORATION  IN  AIR  DURING 
COMPRESSION  IN  A  GAS  TURBINE  ENGINE 

BY;  EARL  QUANDT 

I.  Introduction 

A  water  fog  concept  is  being  considered  for  evaporative  cooling  of  the  air  as  it  is 
compressed  in  a  ship  gas  turbine  engine.  The  following  analysis  is  presented  to  clarify  the 
physics  associated  with  liquid  droplet  evaporation  in  this  situation,  to  understand  the  conditions 
affecting  the  cooling,  and  to  identify  any  further  information  required  to  achieve  such  a  concept. 

The  vaporization  of  small  liquid  drops  in  a  warm  ideal  gas  is  controlled  by  the  outward 
motion  of  the  vapor  and  the  inward  flow  of  heat  to  cause  evaporation.  Following  the  standard 
analysis  of  Spalding,  as  given  in  “Principles  of  Combustion”  by  Kuo*,  it  is  assumed  that  the 
process  is  “quasi  steady.”  This  means  that  the  conditions  far  removed  from  the  drop  are  constant, 
and  that  there  are  no  time  varying  terms  in  the  Eulerian  description  of  the  mass  and  energy  flows. 
For  a  spherically  symmetric  drop  we  have: 


Conservation  of  Mass: 

^  _  0 
dr 

(Total) 

Conservation  of  Mass: 
(Vaporizing  Phase) 

"■f-j,  ( 

Conservation  of  Mass: 

(Ideal  Gas  Phase) 

Conservation  of  Energy: 
(Total) 

dCoT  d 
dr  dr 

It  is  assumed  throughout  that  p.j.  ®  constant,  and  that  a  =  D  =  constant,  Yj  +  Y^  =  L  Y,  « Y^. 

It  will  also  be  necessary  to  assume  that  at  the  droplet  surface  the  vapor  pressure/temperature 
relationship  follows  that  of  the  pure  vaporizing  component. 

II.  Analysis  -  (Quasi-Steady) 

Solution  of  the  above  mass  and  energy  equations  is  possible  by  specifying  boundary 
conditions  at  the  drop  surface  and  at  a  large  distance  from  the  surface.  Then,  these  solutions  can 
be  joined  by  the  vapor  pressure  -  temperature  relationship  at  the  surface.  This  gives  the  “quasi- 


*Kuo,  Kenneth  K.,  “Principles  of  Combustion,”  John  Wiley  &  Sons,  1986. 
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steady”  evaporation  rate,  and  the  “lifetime”  of  a  droplet  of  a  given  size  and  mass.  The  steps  in  the 
standard  analysis  are  given  by  integration  of  the  conservation  equations. 

pvr^  =  const  =  —  =  — i 
47t  47t 

(pvr2)Y  =pr- D^i+C, 
dr 

(pvr2)Y  =  pf-  + 

dr  * 

(pvr2)T  =  pr2  ay  +  C, 

Here  the  similarity  of  the  equations  for  Y  and  T  is  clear,  and  one  only  needs  to  insert  suitable 
boundary  conditions  and  normalizations  to  obtain  identical  solutions.  It  should  also  be  noted  that 
only  the  vaporizing  component  contributes  to  the  radial  mass  flow. 

•  Mass  Concentration 

For  the  water/air  system  of  interest,  we  have 

v  =  l  =  v,Y,.v,Y,=  v,Y,(„'iV^) 

•••P  =  P2 

Thus,  for  constant  pressure,  the  mean  density  in  the  gaseous  phase  is  approximately  constant  as 
long  as  Y,  «Y2,  and  T  does  not  vary  too  much. 

Solving  for  the  several  constants,  we  introduce  the  boundary  conditions  at  the  drop  surface, 
r  =  R, 


and,  since  Y,  +Y,  =  1, 
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(\ 

471  \ 


’'iB*  v,dr 


..J 


o 

=  “^-C 
47:  ' 


Now,  we  recognize  as  proportional  to  the  mass  flux  of  the  ideal  gas  phase,  which  is  zero,  so  C, 
=  m/47r,  and 


_D  ^ 

dr  r  =  R 


=  (I-Y,J 


So,  to  normalize  the  terms  we  let 


b(r).:^ 

‘  ^  IR 

Therefore  we  have  db/dr  I  j.  _ 

(r^pv)b  =  (r^pv)  —  ^ 

V  dr  1-Y.„ 


So  that  at  r  =  R, 


(r^pv)  (b-l)  =  - 


and  in  general 


-..nndb 


(r^pv)  (b-b^+l)  =  r  pZ) 


This  equation  can  be  solved  by  seperating  the  variables  and  integrating  from  r  =  «>  to  r. 


dr_  f^_  db 

^  J  (pr^v)  (b-b^+l 

oo  •'oo 

Yi.I  "1  -  in 

\r  r.  y  (prv) 


-EE.  In  (  b-b,-H 
(pr=v)  Vb.-b,+t 


Therefore 


(b-b,+l) 

(b„-b,+l)-'' 
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and  when  r  =  R,  b  =  b^,  and 

(b.-b^+l)  =  e 

so  that 


b,.b.  =  -e 


Rv„/D 


=  1  -e 


n^47cpD 


and 


b^-b(r)  =  1  -  e 


T^AKpD 


Y.(»-)-Y.r 

1-Y„ 


Thus  we  have  the  mass  concentration  of  Y,  defined  in  terms  of  radial  location  and  evaporation 
rate.  It  is  of  interest  to  note  that  Y,^  is  not  1 .0  at  the  drop  surface. 


Temperature  Distribution 


Turning  to  the  thermal  situation,  the  equation  solution  is  similar,  and  is  related  to  the 
boundary  conditions  at  the  droplet  surface  and  for  large  distances  from  the  drop.  Since  at  the 
surface  the  heat  conduction  through  the  stationary  ideal  gas  phase  provides  the  energy  for 
evaporation,  at  r  =  R,  T  =  T^,  and 


=  R 


-  hvVR  P 


R 


Thus  at  the  surface 

(p  V  r-)  T^  =  p  R^  a  h^  p  v^/?i  +  C, 


Substituting  the  thermal  constant,  C^,  into  the  temperature  differential  equation  gives 

(T  T  +  <1E  _  (xr^  dT 

Cp/  V  dr  Vj^  R^  dr 


Seperating  variables  and  integrating  from  r  =  «>  to  r  =  r  gives; 

ln(T-T„-hh,cJ  -  v,R^  1  -  A 
T.-T.  +  h/c^  -  a  \r  rj 
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Thus  at  r  =  R,  T  =  and 


(T  -T_)  h  /c  +  1  =  e  ''r 
and 

(T  -T_)  h  /c  =  e  - 1  = 

such  that 

(T  -  +  h/c^ )  _  -  m/4;tpar 

T  -T„  +  h  /c  ■  ^ 

oo  R  V  p 

and 

"'■^R  _  e  ^TCap  yR  ry  _  J  _  ^IR  '  ^l(r) 

Therefore,  if  a  =  D,  the  temperature  and  vapor  mass  fraction  solutions  are  equivalent,  if  one 
knows  the  constant  mass  evaporation  rate. 


Vapor  -  Liquid  Eouilibria 


Although  the  equations  for  temperature  and  vapor  mass  fraction  give  similar  functional 
solutions,  we  still  need  an  additional  relationship  between  T  and  Y,  to  give  the  absolute  values. 
This  can  be  obtained  from  the  condition  that  at  the  droplet  surface  the  mass  fraction  of  the  vapor 
is  given  by  its  partial  pressure  for  that  temperature. 


P,.  ",  W,  P„W,  P,(T,.)W, 
"•■p"nW-p,.W"  p,  W 


using  the  Clausius  Clapeyron  relationship. 


it  is  possible  to  express  the  vapor  mass  fraction  at  the  surface  in  terms  of  the  droplet  surface 
temperature.  For  a  reference  condition  on  the  vapor-liquid  equilibrium  line  of  100°C  for  water, 
the  reference  pressure  is  1  atm,  so 


h^.*556cal/g  R  »  0.073  cal/g^k  TR^f=373‘’k 
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and  we  let  T  =  373  +  6t  ,  so  that 

P,  =  e  :07ff373  ('  ■  T^t)  “ '  ““ 

Thus, 

Y-:^ZMe 

'R  ~  W 

so,  we  now  have  three  equations  in  Y,^,  and  or  rn,  so  there  will  be  a  closed  solution  for  R, 
P^,  and  the  properties. 

III.  Application 


In  order  to  illustrate  the  solution  procedure  we  will  examine  the  evaporative  cooling 
behavior  for  a  water  spray  injected  into  the  compressor.  The  initial  conditions  will  be  air  at  100°F 
and  50%  relative  humidity.  We  will  look  first  at  the  evaporation  rate  for  a  single  stage  with 
pressure  ratio  of  1 .3. 

Id 

^^  =  (1.3)  r  *1.08  To  =  311K,T,  =336K&T  =  324K 

At  100°F,  the  vapor  pressure  of  water  is  0.95  psi,  so  that  50%  RH  corresponds 
to  p  =0.475  *  78°F,  and 


Y 


I  «» 


0.475  .  18 
14.7  29 


=  0.02g/g, 


Adjusting  the  Clapeyron  equation  to  100°F  =  31  IK,  p^^^  =  0.95  psi,  h^  =  1040  BTU/LB  =  578 
cal/g 


P,„  =  0.95e 


25.5 


Now,  with  pressurization,  the  water  vapor  initially  in  the  air  is  also  pressurized  to  *  0.63  psi,  or  a 
saturation  temperature  of  87°F.  So,  any  temperature  above  87°F  is  available  to  cause  evaporation. 
At  the  same  time,  the  temperature  rise  of  the  isentropic  compression  has  made  it  possible  for  the 
air  to  hold  0.1 1  gwater/g,  so  the  new  RH  is  18%.  Therefore  T  of  87°  is  a  convenient  reference  for 
the  vapor  pressure,  since  5t=0  results  in  no  concentration  potential  for  Y,. 
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8-L  a  ^ 

<!<  Pl  Vs 

SO  that  the  time  required  to  evaporate  a  drop  size  of  is, 

t  =Px  2a:  f  hJs?i 

®  p  8a  \  AT  / 
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Thus  if  a  =  5  X  10’  mVs  and  =  10’ m,  i.e.  10  microns. 

1.7  X  10^  •10-'»  =  34ms 
°  5x  10 

From  this  baseline  we  may  consider  two  changes,  lower  inlet  humidity  and  behavior  at  the 
higher  pressure  stages  of  the  compressor.  If  the  inlet  humidity  was  zero,  the  6t  will  be  negative  at 
about  -6K,  so  that 

=  0.0 13  and  0.02 1  X  0.47  =  0.0 10 

If5t  =  -4 

_28_  =  0.012  and  0.021  x  0.65  =  0.014 
2312 

So,  a  solution  is, 

6t  =  -5  &  AT/h/c  =0.0125 

V  p 


So  that 

t  =  ^  0^^)  :=  20  ms 

”  p  8a  AT 

Therefore  zero  inlet  humidity  will  act  to  increase  the  evaporation  rate  by  40%. 

In  the  later  stages  of  the  compressor,  where  the  pressure  ratio  might  be  10,  we  have 
^2  =  ( 1 0)°  286  =  2.0  T,  =  622  K 

o 

and  because  of  the  inlet  humidity,  =  4.75  psi,  Tsat  =  160°F  and  h^  =  478  cal/g 
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Therefore,  referring  the  temperature  to  344K,  (160°F) 

p,(6t)  =  4.75e^^^^^^'^^  &  Y,^  (50  =  0.020^^^(144) 

(T„-344)  -  (T„-344)  Y„  -  _  Y„-0.02 

4x478  "  1-Y,„  "  0.98 

Again,  finding  the  solution  to  the  surface  temperature,  surface  vapor  pressure  relationship, 

»  0.02  (0.055  5t)  /0.98 

1912 

But,  if  5t  is  more  than  10°C,  the  small  exponent  approximation  is  not  good,  and  a  more  exact 
interation  is  needed.  Doing  this,  at 

8t  =  35°C  243  =  38.2  (7.5-1)  =  248 

Hence  5t  =  35°C  is  a  solution  that  matches  with  Y,^  at  the  droplet  surface,  and  ^T/h/Cp  = 
0.127. 

Therefore,  for  the  high  pressure  stages  of  the  compressor; 


p,  D  2  (h  /c  )  - 

t  =  _ 0_  '■  V  P''  =r  2 

°  p  8a  ‘^^T 


So,  the  drop  will  evaporate  10  times  faster  in  the  higher  pressure  and  higher  temperature  regions 
of  the  compressor.  It  might  be  noted  that  the  inlet  humidity  will  not  have  a  great  effect  on  this 
result.  Also  a  droplet  temperature  of  380°K  =  225°F,  the  vapor  pressure  on  the  surface  is  19  psia, 
so  Y,^  =  0.08  from  the  steam  tables,  andY,^  =  1.6  x  0.08  from  the  above  Clausius  Clapeyron 
approximation.  Thus,  if  steam  tables  were  used  for  the  iterative  solution,  a  somewhat  higher  5t 
would  result.  This  will  not  have  a  large  effect  on  the  conclusions  regarding  the  evaporation  rate. 


IV.  Closing  Observations 


For  evaporative  cooling  of  the  air  flow  through  a  gas  turbine  compressor  by  water  fog 
injection,  the  following  comments  may  be  made; 


1 .  A  set  of  analytics  is  available  to  compute  an  approximation  to  the  evaporation  cooling  process. 
Solutions  depend  on  a  number  of  assumptions  concerning  simple  transport  properties,  and  quasi¬ 
steady  behavior,  but  are  educational  first  approximations  to  the  true  behavior. 
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2.  The  solutions  show  the  evaporation  rates  to  be  very  droplet  size  and  temperature  sensitive.  For 
standard  gas  turbine  inlet  conditions,  and  10  micron  diameter  drops  at  the  inlet,  it  appears  that, 
for  an  approximate  particle  total  through  transport  time  of  10  ms,  the  cooling  effect  will  be 
delayed  to  the  later  stages  of  compression. 

3.  For  a  100°F  day  with  50%  humidity,  the  mass  fraction  of  water  in  the  inlet  air  is  0.02,  which 
corresponds  to  a  total  water  mass  flow  rate  about  equal  to  the  full  power  fuel  flow  rate.  Adding 
an  equal  amount  of  water  to  the  inlet  as  a  fog  would  provide  a  cooling  capacity  of  40°K.  From  an 
engine  power  output  consideration,  it  would  be  most  beneficial  to  achieve  this  cooling  as  early  in 
the  compression  step  as  possible.  But,  evaporative  cooling  anywhere  should  be  attractive,  the 
only  disadvantage  being  the  large  amounts  of  pure  water  required. 

4.  Water  drops  larger  than  10  microns  may  pass  through  the  compressor  and  into  the  combuster, 
eventually  giving  additional  work  for  additional  fuel.  Cycle  analysis  will  be  needed  to  clarify  the 
rematched  engine  changes  due  to  cooling,  and  to  liquid  carry  over  into  the  combuster. 

5.  With  the  large  amounts  of  water  addition,  existing  data  suggests  that  NOx  reductions  will  be 
quite  significant,  as  already  demonstrated  by  the  WLN  units  currently  used  by  land-based 
utilities. 

6.  It  seems  reasonable  to  explore  further  the  thermodynamic  path  of  droplet  evaporation,  and  gas 
cooling,  in  an  axial  flow  compressor.  Such  parameters  as  transfer  coefficients,  droplet  size 
distribution,  blade  angle  variations  between  wet  and  dry,  gas  temperatures,  etc.  need  to  be 
understood. 

In  summary,  the  concept  of  water  fog  evaporative  cooling  of  the  inlet  air  for  a  ship  gas 
turbine  appears  feasible  and  attractive  from  a  theoretical  standpoint.  The  definition  of  a  practical 
set  of  shipboard  equipment  and  operating  techniques  using  this  concept  should  be  the  next  step 
for  assessment  by  the  Navy. 


10 
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